Probing the non-linear structure of general relativity with black hole binaries 
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Observations of the inspiral of massive binary black holes (BBH) in the Laser Interferometer Space 
Antenna (LISA) and stellar mass binary black holes in the European Gravitational- Wave Observa- 
tory (EGO) offer an unique opportunity to test the non-linear structure of general relativity. For a 
binary composed of two non-spinning black holes, the non-linear general relativistic effects depend 
only on the masses of the constituents. In a recent letter, we explored the possibility of a test to 
determine all the post-Newtonian coefficients in the gravitational wave-phasing. However, mutual 
covariances dilute the effectiveness of such a test. In this paper, we propose a more powerful test 
in which the various post-Newtonian coefficients in the gravitational wave phasing are systemati- 
cally measured by treating three of them as independent parameters and demanding their mutual 
consistency. LISA (EGO) will observe BBH inspirals with a signal-to-noise ratio of more than 1000 
(100) and thereby test the self-consistency of each of the nine post-Newtonian coefficients that have 
so-far been computed, by measuring the lower order coefficients to a relative accuracy of ~ 10~^ 
(respectively, ~ 10~*) and the higher order coefficients to a relative accuracy in the range 10~*-0.1 
(respectively, 10^^-1). 

PACS numbers: 04.25.Nx, 95.55.Ym, 04.30.Db, 97.60.Lf, 04.80.Cc 



Binary pulsar observations provide one of the most 
stringent methods to test the strong field regime of grav- 
ity in general relativity (GR) and its alternatives [J. The 
test is possible since the orbital dynamics of the binary 
is relativistic enough to allow the measurement of ef- 
fects due to gravitational radiation damping at the post- 
Newtonian order (v/c)^ . Binary pulsar measurements are 
performed by fitting the pulse arrival times to a rela- 
tivistic 'timing' model [llljj which is a function of the 
Keplerian parameters (orbital period, eccentricity and 
the projected semi-major axis of the pulsar orbit) and 
post-Keplerian (PK) parameters (the periastron advance, 
time-dilation and secular change of the orbital period). 
Two more PK parameters, related to the Shapiro-delay 
caused by the gravitational field of the companion, can 
be measured if the orbit is seen nearly edge-on. Different 
theories of gravity have different predictions for the val- 
ues of the PK parameters as functions of the individual 
masses of the binary constituents mi and m2. Thus, a 
measurement of three or more PK parameters facilitates 
a test by requiring consistency, within the observational 
errors, in the estimation of the masses of the two bod- 
ies as determined by the various parameters. The most 
rigorous test possible so far is with the most relativis- 
tic binary pulsar PSR J0737-3039 Q- Observed almost 
edge-on, it permitted the measurement of five PK param- 
eters, which together with an additional constraint from 
the measurement of mass-ratio, determine and check the 
consistency of the masses of the two pulsars in the 7711-7712 
plane |3l- 

Although radio binary pulsars are capable of testing 
certain lower post-Newtonian (PN) order general rela- 
tivistic effects, such as the advance of the periastron 
and the quadrupole approximation to the generation of 
gravitational waves, they will, unfortunately, not be able 



to probe the strong field non-linear effects, such as the 
tails of gravitational waves |j| . This is because the PN 
expansion parameter is of order v ^ lO"'^ - far too 
small for the effects that first appear at higher post- 
Newtonian orders to play a significant role in radio ob- 
servations of binary neutron stars. Space- and ground- 
based gravitational wave detectors, such as the Laser In- 
terferometer Space Antenna (LISA), Laser Interferom- 
eter Gravitational- Wave Observatory (LIGO), VIRGO 
and European Gravitational-Wave Observatory (EGO), 
will observe compact binary neutron stars and binary 
black holes (BBH) in the last stages of their non-linear 
evolution, during which the parameter v is two orders 
of magnitude larger (v '^ 0.2-0.4) than it is for current 
radio observations of such systems. For some of the rare 
(about once per year) inspiral events observed by LISA 
(EGO) the amplitude signal-to- noise ratio could be as 
large as 3,000 (100). Such high SNR events will allow us 
to measure the parameters of the signal and the source 
quite accurately, thereby allowing tests that were not fea- 
sible earlier. Different tests of GR have been proposed 
by various authors using GW observations of the inspi- 
ralling compact binaries |a,|3jl3 ^^^ contrasted with the 
binary pulsar observations [g- These tests would neces- 
sitate an accurate parameter extraction J9j scheme using 
the highest PN order waveform available. 

The GW 'phasing formula' is very close in spirit to the 
'timing formula' used in the binary pulsar observations. 
The timing formula, (j)n^^=^T[tn, Pi], connects the rota- 
tional phase (/)„ of a spinning pulsar to the time-of-arrival 
t„ of the radio signal and a set of Keplerian and PK pa- 
rameters pi={p^,p^^}. Similarly, a precise model for 
GWs from a compact binary will need accurate informa- 
tion about the continuous evolution of the GW phase. 
Schematically, the phasing formula reads (j)'^^=Fp[t,qi] 



where, in Einstein's theory, qi carry the information of 
the source via functions of the individual masses and 
spins. The phasing formula consists of different PN pa- 
rameters g,, similar to the PK parameters of the timing 
formula, and is currently available up to relative 3.5PN 
order i.e., Oitf) [lOllllj . In the present paper we propose 
and explore an interesting possibility of testing General 
Relativity with the high-SNR GW observations of BBH 
inspirals by LISA and EGO. The proposed test is similar 
in essence to the binary pulsar test, but in a stronger and 
dynamic regime of gravity. Using the two lowest order 
PN coefficients g,; as basic variables to parametrize the 
waveform and choosing the other PN coefficients as 'test' 
parameters, one at a time, it is possible to perform many 
consistency checks of the PN coefficients in the mi-m2 
plane. In the rest of the paper we investigate this possi- 
bility in greater detail. 

Binary black holes in close orbit around each other are 
highly relativistic and mandate the inclusion of higher 
PN order terms in their description. Gravitational waves 
emitted during the inspiral phase comprise a variety of 
terms arising from the non-linear multipole interactions 
as the radiation propagates from the source to the far- 
zone |l2j . These non-linear interactions lead to the phe- 
nomenon of tails at orders 1.5PN and 2.5PN (propaga- 
tion not only on but inside the light cone as well) and 
tails-of-tails at 3PN. For spinning binaries, there also ex- 
ist effects of spin-orbit and spin-spin coupHngs at 1.5PN 
and 2PN, respectively. These effects are imprinted in 
the emitted gravitational radiation and can be extracted 
by matching the detector data with an expected grav- 
itational waveform, often called an optimal filter or a 
template. The template itself can only be computed us- 
ing post-Newtonian theory in which the various physical 
quantities relevant to the emission of gravitational waves 
are expanded in an asymptotic series in the small param- 
eter V - the characteristic velocity in the system [l^ . An 
important feature of the PN expansion is the presence of 
log-terms i)™(lnw)", where m and n are integers. General 
relativity is incompatible with a simple Taylor expansion 
in only powers of v. For instance, currently, the expan- 
sions of the specific binding energy E and gravitational 
wave flux J-^ are known to order 
and given by 
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(i.e. 3.5PN order) 
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where, rj = miTO2/Af^, is the symmetric mass ratio in 
terms of the total mass M = rrii + 7712 and where the 
coefficients Ek and Tk can be found in Ref. [iJI- Note 
the presence of the log-term at order v^ in the expression 
for the flux. 

To understand how we might test the non-linear struc- 
ture of general relativity let us begin with the Fourier 
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FIG. 1: The signal-to-noise ratio for stellar mass binary black 
holes (BBH) in Advanced LIGO and EGO and supermassive 
BBH in LISA for equal mass binaries at a distance of 200 
Mpc (for EGO and Advanced LIGO) and 2 = 1 (LISA). In 
the case of LISA we assume that the signal is integrated for a 
year (last year before coalescence) and in the case of EGO we 
assume that the signal is integrated over a bandwidth from 10 
Hz until the binary reaches its innermost circular orbit. The 
masses of supermassive BBH in the case of LISA have been 
scaled down by a factor of 10*. 



domain representation if(/) of the signal from a binary 
at a luminosity distance Dl [l^ consisting of black holes 
of masses ra\ and mi: 



H{f) 



CM^I^ 




Dltt^/3 



*(/)+i7r/4 



(3) 



where A4 = rf/^M is the chirp mass and < C < 1 
is a constant that depends on the relative orientation of 
the detector and source with a root-mean-square value 
of 2/5 when averaged over all sky locations and source 
orientations. The phase ^(/) is given by 

7 
*(/) = 2^/t, - $^ + ^ [V,, + V^,, In/] /(''^-5)/3. (4) 
fc=o 

Here tc and <I>c are the epoch of merger and the signal's 
phase at that epoch, respectively. The non-zero coef- 
ficients in the PN expansion of the Fourier phase can 
simply be read off from q;„s in Eq. (3.4) of Ref. j^l- The 
non-vanishing coefficients of the log-terms up to 3.5PN 
areV5/=-if + |ii^ and V'6i=-i^(7rM)i/3. For com- 
pleteness, ijjq=^^{tt MY/^ - ipei log/. We have a total 
of nine post-Newtonian parameters, seven of these are 
the coefficients of w" terms for n = 0, 2, 3,4, 5,6, 7, and 
two are coefficients of v" ln(t!) terms [l^ for n = 5, 6, but 
each of these parameters depends only on the masses of 
the two black holes for the nonspinning case in GR. 

Before proceeding with the description of our work, let 
us summarise the assumptions implicit in the analysis, 
the justification for doing so and their possible implica- 
tions. As in most works on this subject, to demonstrate 
the 'principle' of the proposed method we neglect the 
effects of spins and eccentricity. What will change on in- 
cluding these additional parameters is the accuracy of the 
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FIG. 2: Plot showing the relative errors /S.tpT / '4't , in the test parameters ^t = il^3,tp4,tp5i,il^e,^6i,'<p7, SiS a function of the total 
mass M of a supermassive BBH at a redshift of 2; = 1 observed by LISA (right panel) and of a stellar mass compact binary at 
a distance of Dl = 200 Mpc observed by EGO (left panel). The rest of the details as in Fig0 



test. The spin effects are relevant only when one of the 
black holes is much smaller than the other and/or when 
the black holes have their dimensionless spin angular mo- 
mentum close to unity. It is not clear that astrophysical 
black holes, especially the supermassive ones, will be ex- 
treme Kerr. Except in cases where both BHs are extreme 
Kerr (or close to it) spin effects are less important for the 
proposed tests since we have considered black holes of 
comparable masses in our study. The issue of eccentric- 
ity, especially for certain LISA sources, is a complex issue 
depending on the astrophysical scenario related to forma- 
tion mechanisms of the binary. Our neglect of eccentric- 
ity in these cases is a simplifying assumption at present. 
Finally let us comment on the use of the so-called re- 
stricted PN waveform in this work. Not merely in con- 
nection with tests that have been proposed but more se- 
riously in most works related to the detection problem 
in GW data analysis, the late inspiral and merger part is 
ignored in the first instance. One begins by using state- 
of-the-art restricted PN inspiral templates. Restricted 
PN waveforms will only bring new variety (higher har- 
monics) without increasing the number of parameters; a 
full test should definitely use the full waveform. Includ- 
ing PN amplitude corrections could improve the tests and 
this is what we are doing as a follow-up of the present 
analysis. By the time LISA and EGO operate there could 
be reliable merger waveforms that can be included in the 
phasing and this would make this test more robust. 

Given a high SNR binary black hole event one can, 
in principle, make a model-independent measurement of 
the above PN coefficients by accepting those values that 
best fit the data as our estimates. A procedure in which 
all the parameters 6 = {tc,^c,4>k,4>ki), k = 0,2,..., 7, 
are independently varied to obtain the best possible fit 
of the signal to the data subjects general relativity to 
the most stringent test possible. In a recent paper, we 
explored the power of such a test to determine all the 
known coefficients to a relative accuracy of 100% or bet- 
ter [13. However, this is by no means the most powerful 
test. This is because the covariances between the var- 
ious parameters enhance the errors in their estimation. 



thereby diluting the effectiveness of the test. 

In the present paper we have studied the accuracy with 
which we can measure the PN coefficients by treating at 
a time only three of the nine ipk coefficients to be inde- 
pendent and taking the rest as functions of two of the 
three parameters. Thus, once a high SNR event is iden- 
tified, we suggest to fit the data to a template wherein 
three terms in the PN expansion, rather than just two as 
in detection problem, (or all the PN terms as proposed in 
Ref. [13), are treated as independent parameters. More 
precisely, in Einstein's GR, the tests consist in treat- 
ing the parameters -00 and 1^2 as the fundamental ones 
from which we can measure the masses of the two black 
holes by inverting the relationships i/'o = V'o('Tiij'7i2) 
and i/'2 = '02('7ii, 7712), and asking if the measurement 
of a third parameter, say ipei = ipei{mi,m2), is consis- 
tent with the other two. Instead of the pair {^q, ^2) one 
can, in principle, equally well take any other pair to be 
the fundamental set. The parameters tpo and -02, being 
lower order coefficients, arc best determined as compared 
to the others and constitute our favoured pair |ia |. 

We shall consider the estimation of parameters in the 
ground-based EGO and space-based LISA, using covari- 
ancc matrix, for which we assume the noise PSDs as given 
in Ref. [13 and |23| , respectively. Wc shall take the fun- 
damental parameters to be ipo and -02 in addition to the 
usual extrinsic parameters tc and <I>c. We shall take the 
test parameter ipT to be in turn '03, . . . , 07, ip5i and 06;. It 
should be noted that there is no test corresponding to the 
term involving 05 since it has no frequency dependence 
and simply redefines of the coalescence phase $c. For 
ground-based detectors, Advanced LIGO and EGO, the 
parameters include tc, $c and the three ^j's. For LISA, 
on the other hand, the results correspond to the case of a 
single detector but with amplitude modulation caused by 
the motion of the detector relative to the source. In this 
case our Fourier domain waveform will have amplitude, 
phase and frequency modulations due to the orbital mo- 
tion of LISA and we use the waveform given in Ref [21| . 
Thus for LISA, in addition to the three parameters 
related to our tests we also have the luminosity distance 



and the four angles related to the source's location and 
orientation. 

The power of the tests depends on the SNR achieved 
for the source. In Fig. ^ we have plotted the SNR in 
LISA, EGO and Advanced LIGO lol, for BBH binaries at 
a distance of z = 1 for LISA and a distance of Dl = 200 
Mpc for EGO and Advanced LIGO. In the case of EGO, 
we consider stellar mass BBH of equal masses with the 
total mass in the range IMq to 4001/©, while in the case 
of LISA the mass range is from 10^ Mq to IO'^Mq, but 
scaled down by IC^ so as to fit all the curves in the same 
plot. While the SNR in EGO can reach several lOO's for 
sources that it might observe every once in a year, in the 
case of LISA the SNR could be several 1,000's for the 
supermassive BBH sources that it is expected to observe 
about once per year. The SNR's in both LISA and EGO 
are large enough for the tests to be very powerful probes 
of the PN coefficients and the non-linear effects of GR. 

The lowest order parameters "00 and "02 are measured 
with the smallest errors. In the case of LISA the errors 
for a source at z = 1 are of order 10^^ — 10^^ and in the 
case of EGO the errors for a source at 200 Mpc are of or- 
der lO^"* — lO^''. Fig.|21plots the relative errors Aipx/i'T 
for various parameters tpT as a function of the total mass 
M. From the plots, it is clear that the proposed tests 
can be performed effectively with all 'i/'fe's, especially in 
the case of LISA. This is another reason why LISA is 
such an important mission. All the test parameters, in- 
cluding the log-terms at 2.5PN and 3PN order, can be 
estimated with fractional accuracies better than 10~^ in 
the case of LISA for massive BBH binaries with the total 
mass in the range lO'^-lO^ Mq, and with fractional ac- 
curacies better than 100% in the case of EGO for stellar 
mass BBH binaries with the total mass range 2-lOAf0. 
This demonstrates the exciting possibility of testing the 
non-linear structure of general relativity using the GW 
observations by EGO and LISA. A similar analysis in 
the case of Advanced LIGO for sources with the total 
mass ~ IOMq, shows that all the parameters, except ip^i 
and '06; 1 can be measured to a relative accuracy of 100%. 
Thus, though the 3PN log-term cannot be probed with 
Advanced LIGO, the 2.5PN log-term can be tested lead- 
ing to an interesting possibility in the more immediate 
future. 

With reference to Fig. |21 one may wonder why the er- 
ror in "04 is the largest relative to the other, higher order, 
^'s. We believe that there are several reasons for this 
odd behaviour: recall that the PN terms in the Fourier 
phase are given by ii)kf'^^~^^^^ ■ When k — 5, there is 
no dependence on frequency and when /c = 4 the term 
varies very slowly as f^^^^. Therefore, terms close to 
k = 5 are likely to suffer from large variances since the 
frequency dependence of the corresponding term is weak. 
Although one might expect 06 also to suffer from large 
relative errors, the fact that in this case the term in- 
creases with frequency as f^'^ contributes to making it a 
more important term than -04. We also observe that ip^ 
has significantly larger covarianccs with -00 and and 02 
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FIG. 3: Plot showing the regions in the »7ii-m2 plane that 
correspond to 1-cr uncertainties in the test parameters ■0t ~ 
ip3,i^4,,tp5i,4'6,'tp6i,'>l'7 for a (10^, 10^)Mq supermassive black 
hole binary at a redshift of 2; = 1 as observed for a year by 
LISA. (Note that the 1-a uncertainty in ■03 is smaller than 
the thickness of the line.) 



which adds to its poor determination. 

In Fig 121 we have depicted the power of the proposed 
test in the mi-m2 plane. We present the uncertainty 
contours, with 1-a error bars, associated with the differ- 
ent test parameters in the 7711-7712 plane, when ^0 and 
-02 are used to parametrize the waveform and in the case 
of LISA. The parameter ■06; is much better determined 
by LISA than EGO, as one would expect. This figure is 
an explicit demonstration of the efficacy of the proposed 
test and the accuracy with which the future GW obser- 
vations of BH binaries by EGO and LISA can test GR in 
its strong field regime. 

As mentioned earlier, the spin and angular parameters 
add a lot of structure to the waveform which contain ad- 
ditional information that can be extracted and more tests 
conducted. Covariance between the old and new param- 
eters is likely to increase the error boxes but the tests 
become more demanding as a result of seeking consis- 
tency amongst a greater number of parameters. Future 
studies should look into the more general case incorpo- 
rating the effects of spin and systematic effects of orbital 
eccentricity that could affect the tests, and more interest- 
ingly, go beyond the restricted waveform approximation 
by incorporating the amplitude corrections |22 | to the 
GW phasing. 

We conclude by discussing the extent to which we can 
extend the current proposal to discriminate between dif- 
ferent theories of gravity such as massive graviton theo- 
ries and scalar-tensor theories 0, 123 ■ The limitations of 
GW phasing to quantitatively discriminate between al- 



ternate theories of gravity has been critically discussed 
in |3| and should be kept in mind. For the massive gravi- 
ton theories the IPN phasing term -02 is different and also 
involves the Compton wavelength of the graviton A^ Us- 
ing 4'o and ijj3 as basic variables and 1/^2 as a test ^^ , we 
find that bounds can be set on the value of Xg, modulo 
the neglect of unconiputed higher PN order corrections 
in the theory. Using EGO, which will observe stellar 
mass BH coalescences, we can set a bound on Xg to be 
1.3 X 10"'^^ km whereas with LISA the bounds are as high 
as 7.12 X 10^® km. Scalar-tensor theories like Brans-Dicke 
theory, which predicts dipolar GW emission, have lead- 
ing additional terms in the the phasing formula at a PN 
order lower than in GR. But the dipole GW emission 
is more important for asymmetric binaries than it is for 
equal mass systems. However, for such systems spin ef- 



fects arc also expected to play a crucial role. The present 
paper deals with only non-spinning binaries and we post- 
pone the questions relating to dipolar radiation including 
spin effects to a future work. Once again, these tests will 
be limited by the uncomputed higher order PN contribu- 
tions in the Brans-Dicke theory. 



Acknowledgments 

BRI thanks the University of Wales and Cardiff, U.K. 
for supporting his visit in January 2006 and BSS thanks 
the Raman Research Institute, India for hospitality dur- 
ing August 2005 when part of the work was carried out. 



[1] T. Damour and J.H. Taylor, Phys. Rev. D45, 1840 

(1992). 
[2] T. Damour and N. Deruelle, Ann. Inst. H. Poincare A, 43 

107 (1985); iHd, Ann. Inst H. Poincare A, 44, 263 (1986); 

J. H. Taylor and J.M. Weisberg, Astrophys. J. 253, 908 

(1982). 
[3] A.G. Lyne et al, Science 303, 1153 (2004). 
[4] L. Blanchet and G. Schafer, Class. Quant. Grav. 10, 2699 

(1993). 
[5] L. Blanchet and B.S. Sathyaprakash, Class. Quant. Grav. 

11, 2807 (1994); ibtd, Phys. Rev. Lett. 74, 1067 (1995). 
[6] CM. Will, Phys. Rev. D 50, 6058 (1994); CM. Will, 

Phys. Rev. D 57, 2061 (1998); 
[7] E. Berti, A. Buonanno and CM. Will, Phys. Rev. D 71, 

084025 (2005). 
[8] T. Damour and G. Esposito-Farese, Phys. Rev. D 58, 

042001-1 (1998). 
[9] K.G. Arun, B.R. Iyer, B.S. Sathyaprakash and P.A. Sun- 

dararajan, Phys. Rev. D71, 084008 (2005). 
[10] L. Blanchet, T. Damour, B.R. Iyer, CM. Will and A.G 

Wiseman, Phys. Rev. Lett. 74, 3515 (1995). 
[11] L. Blanchet, B.R. Iyer and B. Joguet, Phys. Rev. D 65, 

064005 (2002); L. Blanchet, G. Faye, B.R. Iyer and B. 

Joguet, Phys. Rev. D 65, 061501(R) (2002); L. Blanchet, 

T. Damour, G. Esposito-Farese and B.R. Iyer, Phys. Rev. 

Lett. 93, 091101 (2004); ibid, Phys. Rev. D 71 124004 

(2005). 
[12] L. Blanchet and T. Damour, Phys. Rev. D 37, 1410 

(1988); L. Blanchet, Class. Quant. Grav. 15, 113 (1998); 

L. Blanchet, Class. Quant. Grav. 15, 89 (1998); 
[13] We use a system of units in which c = G = 1. 
[14] T. Damour, B.R. Iyer, and B.S. Sathyaprakash, Phys. 

Rev. D 63, 044023 (2001); ibid. D 72, 029902E (2005). 
[15] Throughout the paper we assume a cosmological model 

with zero spatial curvature {Qk. = 0, f^A -I- ^m = 



1) and Hubble's constant to be Ho = 70 km s 
Mpc~^. The luminosity distance is given by Dl 
Ijr^ rf — -TPT , where z denotes the redshift 

Hn JO ro_,n_i_-/i3^Q 



11/2 



[nM(i+z . 
of the source. 

[16] In the Fourier domain the log-terms appear at 2.5PN 
order rather than 3PN because of an integration that 
involves 1//. 

[17] K.G. Arun, B.R. Iyer, M.S.S. Qusailah and B.S. 
Sathyaprakash, Class. Quantum Grav. 23, L37 (2006). 

[18] One could alternatively use the entire PN expansion up 
to some k' as a basis, and then check consistency with 
^k'+i (ignoring higher order terms). This would be rel- 
evant for testing alternative theories where deviations 
from GR grow larger with stronger gravity. However, the 
systematic effects from the neglect of known higher order 
PN corrections to GR might contaminate the difference 
between the alternate theory and GR in this case. 

[19] M. Punturo, private communication (fit by A. Sengupta 
and C Van Den Broeck, in preparation). EGO is imag- 
ined to b e a possible third generation gravitati onal wave 
detector. |http://www. astro.cf.ac.uk/geo/ego/| 

[20] L. Barack and C Cutler, Phys. Rev. D 70, 122002 (2004). 

[21] C Cutler, Phys. Rev. D57, 7089 (1998). 

[22] K.G. Arun, L. Blanchet, B.R. Iyer and M.S.S. Qusailah, 
Class. Quant. Grav. 21, 3771 (2004). 

[23] T. Damour and G. Esposito-Farese, Class. Quant. Grav. 
9, 2093 (1992). 

[24] Given the theoretical status of the phasing formula in 
the massive graviton case, this choice we implement is 
the most convenient one at present. In a more complete 
theory -03 could in general depend on the Compton wave- 
length as well and there would be not much rationale be- 
hind the choice proposed. One could then use -00 and tp2 
as basic variables and ■03 as a test, as before. 



